Vortices in condensate mixtures 
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In a condensate made of two different atomic molecular species, Onsager's quantization condition 
implies that around a vortex the velocity field cannot be the same for the two species. We explore 
some simple consequences of this observation. Thus if the two condensates are in slow relative 
translation one over the other, the composite vortices are carried at a velocity that is a fraction 
of the single species velocity. This property is valid for attractive interaction and below a critical 
velocity which corresponds to a saddle-node bifurcation. 
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One remarkable result in condensed matter physics is 
the discovery by OnsagerQ that, in a superfluid, the vor- 
ticity can be present along narrow lines with a quantized 
circulation. Indeed the integral j u.ds of the fluid veloc- 
ity u taken along a circuit enclosing a vortex line must be 
a positive or negative integer multiple of h Planck's 
constant, m mass of the particles making the superfluid. 
This brings a striking analogy between the dynamics of 
the quantum vortices of a superfluid and the Kelvin vor- 
tices of a classical inviscid fluid, the quantization being 
present only to specify the value of the circulation, an 
arbitrary quantity in a classical fluid. This analogy be- 
tween classical inviscid fluids and superfluids is at the 
heart of our understanding of superfluid mechanics, be- 
ginning with the Landau two- fluid theory 2] when there 
is no normal fluid. However things are not so simple, just 
because the quantization of the circulation involves ex- 
plicitely the mass of the particles. Therefore, if there is 
more than one species of particles with different masses, 
it is not obvious at all that classical fluid mechanics re- 
mains the right theory to describe the large scale motion 
of this mixture with quantum vortices. This is because, 
in such a mixture, one does not know a priori which mass 
enters into Onsager's circulation condition. 

This question seems to be irrelevant for superfluid He- 
lium 4, because it has no other stable bosonic isotope, 
and mixing it with any other atomic or molecular liq- 
uid is not possible at temperature low enough to observe 
superfluidity. Mixtures of Helium 4 and 3 (a Fermion) 
can remain liquid, although the spin effects in Helium 3 
make the whole picture quite different, but certainly ex- 
tremely interesting from the present point of view. Below 
we look at a situation that can be, presumably, realised 
in atomic vapors, namely a mixture of two bosonic atoms 
(or eventually molecules) 0, 0, @ . 

We consider the following general problem: given that 
there are two species in the same gas, both condensed, 
what are the dynamical properties of the large scale mo- 
tion of this mixture? This problem, as far we are aware 
of, has not been looked at previously with vortices in- 



cluded ( without vortex look at 0,0)- 

By extrapolating what is known in the case of sin- 
gle species condensate, one can think to many relevant 
questions. The normal modes extending to mixtures the 
Bogoliubov spectrum has been studied by while the 
density profiles in harmonic traps has been investigated 
by [S| . When looking to the fluid motion itself, one of the 
most interesting issues there is the behaviour of vortices. 
We assume that the mixture is at zero temperature and 
that each molecular/atomic species of molecular mass nij 
is described by a macroscopic wave function tyj(r,t), a 
complex valued function of the position r and of time t 
with the discrete index j being cither 1 or 2, to denote 
the species under consideration (one could deal as well 
with more than two species). The equation of evolution 
of the coupled \Pj(r, t) j = 1, 2 is a priori an extension^ 
of the familiar Gross-Pitaevskii equation (G-P later on): 
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The above writing is for two coupled equations, with j = 
1 and j = 2. In the interaction term, the index (j ' + 1) 
is computed mod 2: 1 + 1 = 2, 2 + 1 = 1. Lastly the 
interaction real parameters a,j and g are such that the 
mixture is stable against collapse and against separation 
into two phases, one rich in 1, the other in 2. The stability 
depends on the minimum of the interaction part of the 
energy, the volume integral of 
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The mixture is then stable against collapse if ai and 02 
are positive and if aia 2 > g 2 ■ The linear stability against 
demixing is determined by the Bogoliubov spectrum of 
excitation. We obtain it by seeking the dispersion rela- 
tion between the frequency u) and the wave number k 
of the linear perturbations around homogenous state of 
densities pj = \^j\ 2 respectively. 
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For uncoupled condensates g — 0, we retrieve the Bo- 
goliubov spectrum for each condensate. The condition 
for linear stability against demixing (w real for all k) leads 
to the same criterion a\a 2 > g 2 ■ The coupled equations 
(|T|) are Galilean invariant and one can thus consider the 
flow of both condensates at the same constant velocity 
through Galilean boosts of the wave- functions. More- 
over, notice that for 5 = 0, the uncoupled eq. are 
Galilean invariant separately so that one can consider a 
relative constant flow between each species. For weak 
coupling one can then generalize this property and thus 
consider the relative flow of one species with respect to 
the other one. If the condensates are homogenous, such 
flow remains an exact solution of the equations and the 
model allows an extra "superfluid" property, that is the 
two species can flow one into each other without dissipa- 
tion. For inhomogenous condensates like those contain- 
ing vortices for instance, the interaction between the two 
species generates a friction force and the vortex dynamics 
is affected by the presence of the two species. The goal 
of the present paper is precisely to exhibit such motion 
for simple cases. 

From the coupled equations the vortices bear a dou- 
ble integer index, denoting the numbers of phase winding 
of each wavefunction around the core of the vortex. The 
solution of equations for a vortex (^1,0.2) is of the 
form \& j = e~" lEjt e mj Xj{ r )i where (r, 0) are the polar 
coordinates in the plane perpendicular to the vortex axis, 
and where 

HEj = ajpj + gpj+i- 

The real functions Xj ( r ) are solutions of two coupled or- 
dinary differential equations: 
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where the primes stand for the derivative along the radius 
r (x'j = d-Xj/dr). The asymptotic condition are that, at 
r very large, Xj tends to p®, the uniform density of the 
species j, although \j behaves like r nj at r small. We 
restrict later our study on the two dimensional case al- 
though 3D dynamics should reveal interesting behavior 
(unzipping, Kelvin waves...) and is postponed to further 
work. We assume that the multiply charged vortices (at 
least one of the rij > 1) are unstable and decompose 
into separated single charged vortices as it is the case 
in general for uncoupled condensates pi HtT|. Moreover, 
because of the coupling via the term proportional to g 
in the original equations, the vortices of composite index 
(both rij = ±1) can be stable, with a joint zero at the 



same location. However this is not at all certain in gen- 
eral and such composite vortex can decompose into one 
single charged vortex in each condensate, not located at 
the same position. The interaction between vortices be- 
longing to different species, like vortices (0,1) and (1,0) 
are short ranged, because it depends on the density dis- 
tribution near the vortex core (the interaction between 
vortices of the same species are long ranged, because of 
the velocity field decaying as 1/r at large distances from 
the core) [ly. A reasonable guess is to assume that for 
g negative the vortices of different species attract each 
other (whatever their relative sign), although their inter- 
action is repulsive for g positive. This is based upon the 
fact that the 'interaction energy' is, in a first approxima- 
tion (that is for small g), represented by the integral of 
g(pi — Pi)(p2 ~ P2J1 positive (and repulsive) for g positive 
and negative and attracting for g negative. However, 
because of the Hamiltonian structure of the dynamics, 
this instability is very slow since it manifests through 
radiation coming from the vortex acceleration [11]. Our 
numerics is in complete agreement with this point: we 
observe that the two vortices stand at the same position 
for negative g while they describe a slow outward spiral- 
ing relative motion for positive g. 

We introduce here a convenient dimensionless version 
of the model. Rescaling space and time by the factors 
(mirr^aic^) 1 / 4 /h and (aic^) 1 / 2 /h respectively we obtain 
the following set of two coupled equations: 



--fv^ + ^l^f^ + sl^+Dl 2 ^ (4) 
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with ai = 1/0:2 — y%/mi and fix = 1/02 = 
yoi/ct2- For uncoupled condensates (g — 0), we note 
that the vortex solution ^(r) for each condensate is de- 
termined by a single function /: 



*5(r) = ./ p o/(J^i r ) e ^(e-ftP?t) 
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where Cj = ±1 describes the sign of the circulation of 
each vortex and with / real solution of the eauation|l3|: 



f(r) 



+ (f(r) 2 - l)/(r) = 



Consider now the effect of a flow on a composite vortex, 
when g is negative (that is when this composite vortex 
is stable) . Because of the Galilean invariance of the cou- 
pled equations we have only to study the relative flow 
of one species (say 1) with respect to the other one at 
constant speed (say along direction 1, vi = Uiei). The 
numerics shows that at low speed, the single composite 
vortex splits first into two vortices, a (0, 1) vortex and a 
(1,0) vortex. Then the two vortices move together at a 
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speed that is a fraction (explanation below) of the speed 
of the moving species. At higher speeds, the vortices 
split apart, one being carried by the fluid of the same 
species, the other one remaining immobile, unaffected by 
the velocity of the other species. 

To understand this, the simplest thing is to solve the 
original equations Q by perturbation, assuming the cou- 
pling term to be small as well as the uniform velocity of 
species 1, species 2 being motionless. The zeroth order 
solution is a pair of (1, 0) and (0, 1) vortices, the first 
one being located at i"i(t), the other at r 2 (t). Without 
flow speed and without interactions both vortices remain 
where they are. As speed of species 1 and the interaction 
is turned on, one finds by expansion in powers of a unique 
small parameter, the velocity and the strength of the in- 
teraction, the time dependent solution of the equations. 
The equation of motion for ri (t) and r 2 (t) follows from 
a solvability condition in this expansion. The technical 
details of this derivation will be explained in a coming 
article, but the principle is the same as the one used 
in similar problems, and that is thoroughly exposed in 
the book by Pismen0|. There is a solvability condition 
in the expansion because the ground solution is invari- 
ant under translation of the position of each vortex. At 
the first order, one has to solve a linear equation for the 
perturbation to the basic solution, with an inhomogeous 
term coming both from the external speed of species 1 
and from the interaction. This linear equation has, in 
general, no solution because the homogeneous piece has 
a non trivial kernel. Imposing that the expansion can be 
done amounts to impose the orthogonality (in the sense 
of a scalar product in a function space) of the inhomo- 
geneous term with the kernel of the homogeneous equa- 
tion. This yields a pair of coupled equations for the time 
derivative of the two vortex positions. They read: 
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where ri 2 = r 2 — ri, x corresponds to the vector prod- 
uct, e z is the unit vector perpendicular to the 2D plane 
and g ■ V(\r 12 \) is the potential energy of interaction be- 
tween the two vortices, depending on the norm of r 12 
only : 
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the first one for the relative motion between the vor- 
tices, the second one giving momentum conservation. 
Moreover the rescaled potential U is defined through the 
function / only: 
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The equations of motion keep the Hamiltonian struc- 
ture of the coupled G-P equations: 



dt 



= -e z 
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with H = pieie z • (vi x n) + aaf&tf(ri2). 

However, wave radiations coming from any accelerated 
motion of the vortices have to be added to the dynam- 
ics. To account for this dissipative effects (for the vortex 
motion only, the full set of equations being still Hamilto- 
nian), one needs to estimate the radiative terms coming 
from non uniform vortex motions. Such complicated cal- 
culations have been done for a pair of corotating vortices 
in the G-P model and it shows that the dynamics slowly 
deviates from the Hamiltonian dynamics, with decreasing 
value of the energy [llLll^ |. The stability of the composite 
vortex at zero velocity (stability for negative g only) re- 
lies on this argument. Moreover, this effect disappears for 
the non radiating equilibrium states moving at constant 
speed. They can thus be determined by the Hamiltonian 
dynamics, their stability being determined using this ar- 
gument on the radiative losses. The trajectories of the 
Hamiltonian system follow: 
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where if is a motion constant deduced from the Hamil- 
tonian dynamics 10 • 

An amazing consequence of the equations of motion 
is that there is a possible equilibrium solution at 
constant speed of the two vortices, such that the force 
of interaction is balanced by a kind of Kutta-Joukovsky 
force on each vortex. When this is possible, the joint 
velocity of motion is: 



V(\r 12 \) = / dr(|^ 2 °(r - r 12 ,i)| 2 - ^)(|^(r, t)| 2 - p°) €lpl Vl , 

J expx + e 2 p 2 

^ Using the vortex profile ©, the equation of motion a simple looking result . However this joint drift of the 

simplifies into the set of equations: tw0 vor ti ces cannot happen if the flow speed is too large. 

From the equations (0 EJl , we obtain the following rela- 

<in 2 ff£i£2 dU(r\ 2 ) tion between Vi, ri 2 and the potential if the equilibrium 

"dT = ~ Vl 2^T (eiPl + £2P2)e " X dr 12 solution exists: 
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FIG. 1: Trajectories of the relative vortex motion for different 
velocities a) vi < v m , b) v\ — v m and c) v\ > v m , neglecting 
the radiative effects. This trajectories correspond to constant 
value of the Hamiltonian. 
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Vi = — (eipi +e 2 p2)U {r 12 )e z x . (11) 

From this equation, we deduce first that the separation 
vector r 12 is orthogonal to the imposed velocity v 1 and 
that such a solution can only be found for low enough ve- 
locity. Indeed, the interaction potential has a monotonic 
behavior with zero derivative both at zero distance and 
at infinity, so that the derivative has a maximum in be- 
tween which determines as well the maximum value v m 
of the velocity defined in eq. I|llfl . Thus two solutions 
exist for velocity smaller than this maximum. The two 
solutions collapse for the maximum velocity and no more 
equilibrium positions exist above. Figure shows the 
Hamiltonian structure of the dynamics for the four cases 
V\ = 0, V\ < v m , Vi = v m and vi > v m . The stabil- 
ity of the equilibrium solutions depends then also on the 
sign of g. For g > every solution is linearly unstable 
and no steady flow can be obtained. For g < 0, we have 
a saddle-node bifurcation structure: at low speed when 
we have two equilibrium positions, one of the solution is 
linearly stable (the one of smaller ri 2 ) and the other one 
unstable. For v\ = 0, the trajectories are circular and de- 
pending on the sign of g the vortex dynamics corresponds 
to a slow converging (diverging) spiraling motion towards 
(outward) the origin. Notice also this surprising effect: 
consider e 2 — — e\ and p 2 > p\, then the constant speed of 
the two paired vortices is in the opposite direction of the 
imposed flow velocity! Such counterflow vortex dynam- 
ics has been observed in the numerics as well. Moreover, 



in the case e 2 = — ei and p® = Pi the only equilibrium 
solution is for Vi = 0. The experimental consequences of 
this arc simple to state in atomic vapors, where one can 
manipulate the condensate by optical methods in par- 
ticular. It would be interesting also to see if, in super- 
fluid mixtures of Helium 4 and 3, there are two kinds of 
vortices, but it is far less easy, but perhaps not impossi- 
ble, there to move independently the two kind of atoms. 
Another last remark is about the rotating two species 
condensate. With a single species, the equilibrium state 
is a triangular lattice of like sign vortices with a mesh 

1/2 

size of order (t^W) , ^ being the angular frequency. 
Therefore, if the two coupled condensate are subject to 
the same angular speed tt, each bosonic species should in 
the absence of coupling g = exhibit a lattice of vortices 

/ \ !/ 2 

of mesh ( -^-^ J . A non zero coupling induces most 

probably deformations of the two lattices. Even in the 
weak coupling limit, no formal theory seems available for 
this kind of situation. Another instance where coupled 
vortices would be present is the one of nonlinear optical 
fields in the classical approximation^). There, the role 
of time in the G-P equation is played by the direction par- 
allel to the beam propagation and the Laplacian account 
for the 2D variation perpendicular to it. The equations 
of propagation of two parallel light beams of different fre- 
quency in the same nonlinear material are very similar 
to the coupled G-P equations. Therefore, we expect in 
this case the occurence of phenomena roughly similar to 
the one descibed here. 
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